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Abstract 

In this paper, the dynamics of the Chebyshev-Halley family is studied on quadratic polynomials. A singular set, that 
we call cat set, appears in the parameter space associated to the family. This cat set has interesting similarities with the 
Mandelbrot set. The parameters space has allowed us to find different elements of the family such that can not converge 
to any root of the polynomial, since periodic orbits and attractive strange fixed points appear in the dynamical plane of the 
corresponding method. 
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I. Introduction 

The application of iterative methods for solving nonlinear equations f(z) = 0, with / : C — > C, give rise to rational 
functions whose dynamics are not well-known. The simplest model is obtained when f(z) is a quadratic polynomial 
and the iterative process is Newton's method. The study on the dynamics of Newton's method has been extended to other 
point-to-point iterative methods used for solving nonlinear equations, with convergence order up to three (see, for example 

II, Jl and, more recently, [3] and 0]). 

The most of the well-known point-to-point cubically convergent methods belong to the one-parameter family, called 
Chebyshev-Halley family, 



Zn+l 

where 



1 , 1 L f( z n) \ f(Zn) m 

+ 21 -aL f (z n ) J f'(z n y K> 

L f ( Z ) = lMI^l, (2) 
(f(z)f 

and a is a complex parameter. This family includes Chebyshev's method for a = 0, Halley's scheme for a = |, super- 
Halley's method for a = 1 and Newton's method when a tends to ±oo. As far as we know, this family was already 
studied by Werner in 1981 (see [5]), and can also be found in ]6|] and |7|]. Moreover, a geometrical construction in studied 
in HI]. It is interesting to note that any iterative process given by the expression: 

z n +i = z n - H (L f (z n )) , (3) 

where function H satisfies H(0) = 0, H'(0) = | and |if"(0)| < oo, generates an order three iterative method (see Q9I1). 
The family of Chebyshev-Halley has been widely analyzed under different points of view. For example, in flioll . ill ill 



and 1 12], the authors studied the conditions under the global and semilocal convergence of this family in Banach spaces 
are hold. Also, the semilocal convergence of this family in the complex plane is presented in [13]. 

Many authors have introduced different variants of these family, in order to increase its applicability and its order 
of convergence. For instance, Osada in 11411 showed a variant able to find the multiple roots of analytic functions and a 
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procedure to obtain simultaneously all the roots of a polynomial. On the other hand, in JTll . Hill and Jnll the authors got 
multipoint variants of the mentioned family with sixth order of convergence. Another trend of research about this family 
have been to avoid the use of second derivatives (see |1_8] and fl^P or to design secant-type variants (see i2(ilp . 



From the numerical point of view, the dynamical behavior of the rational function associated to an iterative method 
give us important information about its stability and reliability. In this terms, Varona in |21] described the dynamical 
behavior of several well-known iterative methods. More recently, in and |22], the authors study the dynamics of 
different iterative families. 

1.1. Basic concepts 

The fixed point operator corresponding to the family of Chebyshev-Halley described in (fl} is: 

In this work, we study the dynamics of this operator when it is applied on quadratic polynomials. It is known that the 
roots of a polynomial can be transformed by an affine map with no qualitative changes on the dynamics of family ([T|i (see 



023IP . So, we can use the quadratic polynomial p (z) = z 2 + c. For p(z), the operator dU corresponds to the rational 



function: 



_ , . _ z A (-3 + 2a) + 6cz 2 + c 2 (1 - 2a) 



Az (z 2 (-2 + a) + ac) 

depending on the parameters a and c. 

P. Blanchard, in 12411 . by considering the conjugacy map 

h(z)= , (6) 

Z + ty/C 

with the following properties: 

i) h (oo) = 1, ii) h (i^/c) = 0, iii) h (—i^fcj = oo, 

proved that, for quadratic polynomials, the Newton's operator is always conjugate to the rational map z 2 . In an analogous 
way, it is easy to prove, by using the same conjugacy map, that the operator G p (z) is conjugated to the operator O p (z) 

O p (z) = (ho Gp o h- 1 ) (z) = z 3 *~*£~^ z - (7) 

In addition, the parameter c has been obviated in O p (z). 

In this work, we study the general convergence of methods (fTJ for quadratic polynomials. To be more precise (see 
1 25] and 1 26]), a given method is generally convergent if the scheme converges to a root for almost every starting point 
and for almost every polynomial of a given degree. 

7.7.7. Dynamical concepts 

Now, let us recall some basic concepts on complex dynamics (see B27I0 . Given a rational function 7? : C — > C, where 
C is the Riemann sphere, the orbit of a point zq 6 C is defined as: 



z , R{z ),R 2 (z ),...,7T (z ),... 

We are interested in the study of the asymptotic behavior of the orbits depending on the initial condition zq , that is, we 
are going to analyze the phase plane of the map 7? defined by the different iterative methods. 

To obtain these phase spaces, the first of all is to classify the starting points from the asymptotic behavior of their 
orbits. 

A zo £ C is called a fixed point if it satisfies: 7? (zo) = Zq. A periodic point zq of period p > 1 is a point such that 
R p (zq) = zq and R k (zq) ^ zq, k < p. A pre-periodic point is a point zq that is not periodic but there exists a k > 
such that R k (zq) is periodic. A critical point zq is a point where the derivative of rational function vanishes, 7?' (zo) = 0. 

On the other hand, a fixed point zq is called attractor if |7i'(zo)| < 1, superattractor if |7i'(zo)| = 0, repulsor if 
|7?'(^o)| > 1 and parabolic if \R' (zq)\ = 1. 

77ie basin of attraction of an attractor a is defined as the set of pre-images of any order: 

A (a) = {z e C : 7?" (z ) ->a, ?i->oo}. 
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The set of points z E C such that their families {R n (z)} n£N are normal in some neighborhood U (z) , is the Fatou 
set, J- (R) , that is, the Fatou set is composed by the set of points whose orbits tend to an attractor (fixed point, periodic 
orbit or infinity). Its complement in C is the Julia set, J (R) ; therefore, the Julia set includes all repelling fixed points, 
periodic orbits and their pre-images. That means that the basin of attraction of any fixed point belongs to the Fatou set. 
On the contrary, the boundaries of the basins of attraction belong to the Julia set. 

The invariant Julia set for Newton's method is the unit circle S* 1 and the Fatou set is defined by the two basins of 
atraction of the superattractor fixed points: and oo. On the other hand, the Julia set for Chebyshev's method applied to 
quadratic polynomials is more complicated than for Newton's method and it has been studied in [28]. These methods are 
two elements of the family (Q~|i. In the following sections, we look for the Julia and Fatou sets for the rest of the elements 
of the mentioned family. 

The rest of the paper is organized as follows: in Section 2 and 3 we study the fixed and critical points, respectively, 
of the operator O p (z). The dynamical behavior of the family (Q]) is analyzed in Section 4. We finish the work with some 
remarks and conclusions. 



2. Study of the fixed points 

We are going to study the dynamics of the operator O p (z) in function of the parameter a. In this section, we calculate 
the fixed points of O p (z) and in the next one, its critical points. As we will see, the number and the stability of the fixed 
and critical points depend on the parameter a. 

The fixed points of O p (z) are the roots of the equation O p (z) = z, that is, z = 0, z = 1 and 

-3 + 2a ± V5 - 12a + 4a 2 
Z = 2 ' (8) 

which are the two roots of z 2 + (3 — 2a)z + 1 = denoted by s% and S2- 

The number of the finite fixed points depends on a. Moreover, s% = — , so that, these points are equal only if 

Si = S2 = ±1; this happens when 5 — 12a + 4a 2 = 0, i.e., for a = i and a = I. 

For a = 5, si = s 2 = — 1, so z = —1 is a fixed point with double multiplicity. For a = |, si = s 2 = 1, so z = 1, 
has multiplicity 3. 

Summarizing: 

• If a ^ h and a ^ |, there are five different fixed points with multiplicity 1. 

• If a = \, there are four different fixed points: z = 0, z = 00 and z — 1 with multiplicity 1 and z = — 1 with 
multiplicity 2. 

• If a = I, there are 3 different fixed points: z = and z = 00 with multiplicity 1 and z — 1 with multiplicity 3. 

As we will see in the next section, the multiplicity of the fixed points implies different dynamical behaviors. 
In order to study the stability of the fixed points, we calculate the first derivative of O p (z), 

, , s 9 3 (1 - a) + 2z (3 - 4a + 2a 2 ) + 3z 2 (1 - a) 

O' (z) = 2z 2 — i ^ i . (9) 

(l-2(a-l)z) 2 

From we obtain that the origin and 00 are always superattractive fixed points, but the stability of the other fixed 
points changes depending on the values of the parameter a. These points are called strange fixed points. 
The operator O p (z) in z = 1 gives 



2 + a)(2a-3) 




4a -8 


(2a - 3) 2 




2a -3 



(10) 



If we analyze this function, we obtain an horizontal asymptote in \0' (1) = 2, when a — > ±00, and a vertical 
asymptote in a = I. In the following result we present the stability of the fixed point z = 1. 

Proposition 1. The fixed point z — 1 satisfies the following statements : 

i) If J a — ^ I < |, then z = 1 is an attractor and, in particular, it is a superattractor for a = 2. 

ii) If a — 4r \ — h, then z — 1 is a parabolic point. 
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Hi) If \a — ^ > 5, f/zen z = 1 is a repulsive fixed point. 



Proof. From equation ( [Tol l. 



< 1 => 4 



Let a = a-f 

and 
So, 

By simplifying 
that is, 
Therefore, 



2a- 3 

ib be an arbitrary complex number. Then, 



a I < 1 2a - 31 



|-2 + a| 2 = {-2 + a) 2 + b 2 

|2a-3| 2 = (2a- 3) 2 + 46 2 . 
16 (4 - 4a + a 2 + 6 2 ) < 4a 2 - 12a + 9 



46 2 



55 - 52a + 12a 2 + 126^ = 12 ( a - — 

6 



12b 2 



13 



b 2 < 



I Op (1)| < 1 if and only if 



13 



1 

< -. 

- 3 



Finally, if a satisfies |a — i 3 -] > i, then |O p (1)| > 1 and z = 1 is a repulsive point. ■ 

The stability of the other strange fixed points z = Si, i = 1, 2 also depends on parameter a. 

|0£(«<)| = |6-2a|. 

We can establish the following result: 

Proposition 2. The fixed points z = Si, i = 1,2 satisfy the following statements: 

i) If \a — 3| < i, f/ien Si anc/ S2 sfs ftvo different attractive fixed points. In particular, for a = 3, Si ant/ S2 are 
SMperaffracfori. 

n'J T/ 1 |q — 3| = |, f/ze« Si S2 fl^e parabolic points. In particular, for a = I, Si = S2 = 1. 

/n'J /f |a — 3| > 4, f/zen Si and S2 are repulsive fixed points. 

In the following bifurcation diagram (FigureHJ we represent the behavior of the fixed point for real values of parameter 

a. The point z = 00 is not represented. Let us observe that the stability of the fixed points is represented by the thickness 

of the lines: if it is attractive, the line corresponding to the value of this strange point is thicker. So, it can be noticed that 

z — is always an attractor, meanwhile z = 1 is attractive when ^J- < a < | and Sj, i = 1, 2 are attractors only when 

- < a < - 
2 ^ ^ 2' 

3. Study of the critical points 

Let us remember that the critical points of O p (z) are the roots of 0' p (z) = 0, that is, z = 0, z = 00, and 

3 - 4a + 2a 2 ± V-6a + 19a 2 - 16a 3 + 4a 4 



3 (a - 1) 



(11) 



which are denoted by c\ and C2. 



It is easy to prove that c\ = — . Therefore, both critical points coincide only when c\ — C2 = ±1, that is, when 

C2 

-6a + 19a 2 - 16a 3 + 4a 4 = 0. 

1 3 

The roots of this equation are 0, — , — and 2. 

It is known that there is at least one critical point associated with each invariant Fatou component. As z — and 
z — co are both superattractive fixed points of O v (z), they also are critical points and give rise to their respective Fatou 
components. For the other critical points, we can establish the following remarks: 
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Figure 1 : Bifurcation diagram of fixed points 



a) If a — 0, then c\ — c 2 = —1, and it is a pre-image of the fixed point z = 1: O v {— 1) = 1. As z = lis repulsive, 
z = — 1 e J {O p ). So, O p {z) has precisely two invariant Fatou components, A (0) and ^4 (oo). 

b) If a = h, then ci = C2 = — 1 = s% = s% are repulsive fixed points and belong to Julia set. 

c) If a = I, then c\ = c 2 = 1 is a repulsive fixed point and belong to Julia set. 

d) If a — 2, c\ = c 2 = 1. In this case z = 1 is a superattractor, which gives rise to a Fatou component. 

e) For any other value of a E C, there are four different critical points, and we will study their behavior in Section 4. 

In Figure [2] we represent the behavior of the strange fixed points and critical points for real values of a between 1 
and 4. We observe that the critical points Cj, i = 1, 2 are inside the basin of attraction of z = 1 when it is attractive 
(-77- < a < |) and coincide with z = 1 for a = 2. Then, they move to the basins of attraction of s± and S2 when 
these fixed points become attractive (| < a < ^), critical and fixed points coincide for a = 3 and s% and S2 become 
superattractors. 



4 

3 
2 
1 


ft • a 




2.0 2.5 3.0 3.5 4.0 



Figure 2: Dynamical behavior of strange fixed points and critical points for 1 < a < 4 

Moreover, we can see that when a 1, C\ tends to and c 2 tends to 00. This fact explains that O p (z) = z 4 when 
a = 1 (super-Halley's method), and the only superattractive fixed points were and 00. 
Finally, if a —> ±00, C\ tends to and c 2 tends to ±00 and O p (z) = z 2 . 

4. The parameter space 

It is easy to see that the dynamical behavior of operator O p {z) depends on the values of the parameter a. In Figure 
[3j we can see the parameter space associated to family (Q]): each point of the parameter plane is associated to a complex 
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Figure 3: Parameter plane 



value of a, i.e., to an element of family (|T). Every value of a belonging to the same connected component of the parameter 
space give rise to subsets of schemes of family ([TJ with similar dynamical behavior. 

In this parameter space we observe a black figure (let us to call it the cat set), with a certain similarity with the known 
Mandelbrot set (see 1291): for values of a outside this cat set we will see, numerically, that the Julia set is disconnected. 
The two disks in the main body of the cat set correspond to the a values for those the fixed points z = 1 (the head) and 
z = s\ and z — s 2 (the body) become attractive. Let us observe that the head and the body are surrounded by bulbs, of 
different sizes, that yield to the appearance of attractive cycles of different periods. 

We also observe a closed curve that passes through the cat's neck, we call it the necklace. As we will prove in the 
following, the dynamical planes for values of a inside this curve are topologically equivalent to disks. 

4.1. The head of the cat set 

The head of the cat corresponds to the values of parameter a for which the fixed point z = 1 become attractive, that 
is, the values of a such that la — 4r| < i. 

16 13 

In this case, the fixed point z = 1 is an attractor (PropositionQJi and the other two fixed points si and s 2 are repulsors 
(Proposition |2). Depending on the values of the parameter the critical points ci and c 2 have different behaviors around 
z = 1. In Figure|2]we observe the behavior of the critical points c\ and c 2 for real values of parameter a in the interval 
(-^-, 2). For iJ- < a < 2 both critical points are complex. When a = 2, both critical points coincide with the fixed point 
z = 1, so that, it is an superattractor. 

4.1.1. The boundary of the head 

As it has been established in PropositionQ] the boundary of the previous set (la — I = |) is the loci of bifurcation 
of the fixed point z = 1. This fixed point is parabolic on this boundary and yields to the appearance of attractive cycles, 
as it happens in Mandelbrot set when we move into the bulbs (see f29lD . 

In this region, a = ^ + ^e l9 and the operator (|7]) can be expressed as: 




(12) 



The first derivative is: 




2 -63 + 134z- 18e 



- 63z 2 + 56ze w - 18z V s + 8ze [ 



= z 



(-3 + 7z + 2ze i6 f 



and it is easy to check that z = 1 is an parabolic point for all these values of a, since 




2e ie + 1 
2 + e ie 



and \<y (1)1 = 1. 
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Therefore, for different values of 9, we find the different bulbs with attractive cycles surrounding "the head of the cat" 
In this boundary there exist two points that are specially interesting: they correspond to the intersection with the real 
axe: 8 — n and 9 = 0. In the first case, a = the operator has the expression 

O p (z) = z A 



3-5z 



Forthis value of a, the two strange fixed points s,; i = 1, 2 are repulsive. Thepointz = 1 is parabolic (since O p (1) = — 1), 
and it is in the common boundary of two parabolic regions : the elements of the orbit corresponding to an initial estimation 
in one of these parabolic regions, go alternatively from one region to the other while approaching to the parabolic point 
z = 1. If a is close to but lower than z = 1 is a repulsive fixed point (see PropositionQ]). As it happens in Mandelbrot 
set, when we take a value of a in this bulb, an attractive cycle of period 2 appears. 



For the case = 0, a = | and 



O p (z) 



z 



3 



1 - 3z 



Now, the three strange fixed points are the same, z = 8\ = s 2 = 1, and it is parabolic, \0' p (1)| = 1. We know, by the 
Flower Theorem of Latou (see OOfl . for example), that this parabolic point is in the common boundary of two attractive 
regions. The orbits of initial estimations inside each region approach to z = 1 without leaving its region. These attractive 
areas contain the respective critical points c\ = | (ip — VlO) and c 2 = § (3- + vTOj. 

4.2. The body of the cat set 

As we have said before, the body of the cat set corresponds to values of the parameter such that \a — 3| < 5 - In this 

case, 

.3 * ~ 2 (<* ~ 1) 
" ) z ' 



°^) = * \ 2 fa"l 
1 — 1 [a — \) . 



the fixed point z = 1 is a repulsor (PropositionQ]) and s,;, i = 1, 2 are attractors (Proposition|2]l. So, they have their own 
basins of attraction with a critical point in each one (see Figure|2]i. 
We know that, for a = 3, 

z = 1 is a repulsor and Sj = Cj, i = 1, 2 are superattractors. 

4.2.1. The boundary of the body 

Similarly to what happen in the Mandelbrot set, the boundary of the cat set is exactly the bifurcation locus of the family 
of Chebyshev-Halley operators acting on quadratic polynomials; that is, the set of parameters for which the dynamics 
changes abruptly under small changes of a. 

This boundary correspond to values \a — 3| = h, that is, a = 3 + \e . 

The strange fixed points Sj, i = 1, 2 are parabolic. So, the different values of the argument 9 give the bifurcation 
points for the different bulbs surrounding the body of the cat set. 

Two interesting values of a are the intersection between this boundary and the real axe: they correspond to 9 = and 
9 = 7T. If 9 = 0, a = I and 

In this case, z — 1 is a repulsor and the two strange fixed points s\ and s 2 are parabolic. Each of these points is in the 
common boundary of two parabolic regions, the iterations of the orbit of an initial estimation inside one of these regions 
go alternatively from one area to the other, while approaching to the parabolic point. 

In the parameter space, the bulb corresponding to values of a bigger and close to | is the loci of the cycles of period 

2. 

When 9 = n, a = | and this value of the parameter (see Figure [3]) corresponds to the intersection between the 

5 



boundaries of the body and the head of the cat. So, the point a = § is a bifurcation point and the dynamics changes when 

5 
2 



the parameter a varies in a small interval around I and it has been studied in the previous section. 
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4.3. Inside the necklace 

For values of the parameter inside the necklace, by applying Propositions[T]and|2j the only superattractive fixed points 
are and oo; the Julia set is connected but we see that for different values of a the fixed points and oo have one connected 
component in each basin of attraction. As we will see in Proposition [3] if | < a < §, these connected components are 
disks in Riemann sphere. In other cases, they are topologically equivalent to disks. In the following section, we prove 
this statement for \a — 1| < h, although it can be checked numerically that it is also true in the rest of the area inside the 
necklace. 

4.3.1. The region \a — 1| < ^ 
In this area, 

O (z)-z 3 z ~ 2 ^- 1 ) 
u P {z)-z 1 _ 2(a _ 1)2 = 

where a = 1 + re l ° , r < i. 

Proposition 3. If \a — 1| < |, then the dynamical plane is the same as the one of z n . 

Proof. If a G R, the map m{z) — iZ^a-ijz * s a Moebius map that sends the unit disk in the unit disk. There are 
only two basins of attraction, of and oo. We are going to prove that the critical point c\ is inside the basin of attraction 
of and C2 is inside the basin of attraction of oo, that is, we need to prove that 



|ci| 

or equivalently, 



3 - 4a + 2a 2 - \/-6a + 19a 2 - 16a 3 + 4a 4 



< 1, 



3 - 4a + 2a 2 - y/-6a + 19a 2 - 16a 3 + 4a 4 



3 (a - 1) 

' <3|(a-l)| 



For A < a < f , it is easy to prove that 3 — 4a + 2a 2 — yj— 6a + 19a 2 — 16a 3 + 4a 4 > 0. So, we need to consider 
only two cases: 



i) if a > 1, we want to demonstrate that 3 — 4a + 2a 2 — %/ —6a + 19a 2 — 16a 3 + 4a 4 < 3 (a — 1) is verified. By 
simplifying, it is equivalent to (2a — 3) (a — 2) < a (2a — 1) and this is equivalent to 6 < 6a. 



ii) if a < 1, then we need to prove that 3 — 4a + 2a 2 — V— 6a + 19a 2 — 16a 3 + 4a 4 < 3 (—a + 1) and, in an 
analogous way as before, it is easy to prove that 6a < 6. 

Moreover, as I C2 1 = -, — r th en l c 2| > 1- So, c\ is in the basin of and C2 in the basin of oo. The Julia set is the unit 

circle that divides these two basins. 

If a is a complex number, the map m (z) = is not a Moebius map, but it is holomorphic. Let us analyze 

the mapping of unit circle by m. The pole of this map is z* — — and, in this case, \z*\ > 1. 

2 (a — 1) 

2 1 

Let z = x + iy be a complex number such that \z\ = 1 and a = a + ib, then, (a — 1) + b 2 < — . Let us see the value 



of |m (z)\ 



z-2(q-l) 
l-2(a-l)z 



|z-2(a-l)| 2 = (x - 2a + 2) 2 + (y - 2b) 2 



1 + 4 ((a - i f + b 2 ^j - 4ax + 4x - Ayb 



|l-2(a-l)z| 2 = (1 - 2ax + 2x + 2ybf + {-2ay + 2y - 2bxf 
= 1 + 4 ((a - l) 2 + 6 2 ) - 4ax + 4x + 4yb 

We observe that both are equal if and only if b — 0, that is, in the real case. If yb > then |m (z)\ < 1 and yb < 
implies \m (z)\ > 1. As m (z) is a holomorphic function, the image of the unit circle is a closed curve that separates the 
images of the points inside the unit circle from those that are outside it. So, the dynamical plane for the values of the 
parameter inside this range consists on two regions of attraction, A (0) and A (oo) separated by this closed curve. As 
before, by continuity each critical point is in one of these regions. ■ 

Therefore, it follows that the dynamical plane of this operator is, for the given values of a, equivalent to the one of z n . 
In particular, we observe that, for quadratic polynomials: 
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• for a = 1, O p (z) = z 4 , 

• if a = | then O p (z) = z 3 

• and, for a = |, O p (z) — — z 3 . 

Let us remark that dynamical plane associated to every iterative algorithm whose value of a is inside the necklace, is 
topologically equivalent to the previous one. 

4.4. On the necklace 

We focus our attention on the real values of a included in the necklace. In particular, in those which belong to the 
antennas of the cat set. If < a < — and a E K, we move in the left antenna of the necklace in the parameter space, that 

is in the boundary of the cat set. We prove in the following result that, in the dynamical plane associated to the iterative 
methods defined by these values of a, there are infinite connected components of the basins of attraction, corresponding 
to the immediate bassins of attraction and their pre-images. 

Proposition 4. The dynamical plane for values ofaEM. and < a < | consists in two basins of attraction, A(0) and 
A(oo) with infinity pre-images. 

Proof. If < a < — and a E K, we move in the left antenna on the necklace, that is in the boundary of the cat set. 
For these values, the strange fixed points are repulsive, so they belong to the Julia set. Moreover, the critical points verify: 

(3-4a + 2a 2 ) 2 + (j-a (a - 2) (2a - 1) (2a - 3)) ^ 

|Cl|2 = V q / U 2 L = l (13) 

9 (a — 1) 

and, as c 2 = — , Icgl = 1. So both critical points are on unit circle. 

Cl 

Moreover, the operator (0 has a pole in z* — — such that \z*\ < 1; so, there is an image of oo inside the unit 

2 (a- 1 ) i—i 
disk and, by symmetry, there is an image of zero outside the unit disk. So, by the Theorem of Fatou (see B30IO . they have 

infinity basins of pre-images. ■ 

We observe the same dynamical behavior for values of the parameter in the right antenna of the cat. The reason is that 
3 

for — < a < 2 (whose values include the ones of the right antenna) and a E K we can use the relationship ([T3T l. 

The case of a = has the same operator on quadratic polynomial that the ones studied by the authors in 113 111 for other 
different iterative methods. The dynamical plane is similar to the previously described. 

4.5. Outside the cat set 

The cat set, as the Mandelbrot set, could also be defined as the connectedness locus of the family of rational functions 
of Chebyshev-Halley methods. That is, it is the subset of the complex plane consisting of those parameters for which the 
Julia set of the corresponding dynamical plane is connected. All the dynamical behaviors we have studied for values of 
the parameter outside the cat set show disconnected Julia sets. 



5. Conclusions 

We have studied the dynamics of the Chebyshev-Halley family when it is applied on quadratic polynomials. We have 
obtained the fixed and critical points and their dynamical behavior, and we have showed that strange fixed points appear 
which are attractive for some values of a (Propositions[T]and|2]l. This means that, for these values, these iterative methods 
have basins of attraction different of the roots of the polynomial. So, the initial point must be chosen carefully. 

From the parameter plane obtained, we have observed the cat set, with some similarities with Mandelbrot set: the head 
and the body of the cat are surrounded by bulbs. For values of the parameters inside the bulbs, different attractive cycles 
appear. We have also studied the dynamical behavior of the family for values of a inside the necklace and we have shown 
that it is analogous to the dynamical behavior of Newton's method (Proposition 3). For values of a on the antennas of the 
cat the dynamical plane has only two basins of attraction, but these basins have infinitely many components (Proposition 
4). Finally, we have obtained numerically that the Julia set is disconnected for values of a outside the cat set. 
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The cat set is a fascinating creature of complex dynamics. Similarly to the Mandelbrot set, there is a lot of different 
dynamics for this cat. We study some of its properties in this paper, but we are aware that there are plenty of unresolved 
issues, for example, it is connected the cat set? We conjecture that the cat set is connected. 

We have also observed little cats in the necklace. What's about the dynamics for these values of the parameter? Even 
more, where are these cats exactly?. What happens in the antennas for non real values of a? 
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